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ABSTRACT 

We report the results of Suzaku broad band X-ray observations of the two narrow-Hne radio galaxies 
(NLRGs), 3C 403 and IC 5063. Combined with the Swift/B AT spectra averaged for 58 months, we 
are able to accurately constrain their spectral properties over the 0.5-200 keV band. The spectra of 
both nuclei are well represented with an absorbed cut-off power law, an absorbed reflection component 
from cold matter with an iron-K emission line, and an unabsorbcd soft component, which gives a firm 
upper limit for the scattered emission. The reflection strength normalized to the averaged BAT flux is 
i? = r2/27r ~ 0.6 in both targets, implying that their tori have a sufficiently large solid angle to produce 
the reprocessed emission. A numerical torus model with an opening angle of ~ 50° well reproduces 
the observed spectra. Wc discuss the possibility that the amount of the normal gas responsible for 
Thomson scattering is systematically smaller in radio galaxies compared with Seyfert galaxies. 
Subject headings: galaxies: active - galaxies: individual (3C 403, IC 5063) - X-rays: galaxies 



1. INTRODUCTION 

Modern astronomy has revealed that the mass of 
a super-massive black holes (SMBHs) in a galactic 
center has tight correlation with the mass of the 
bulge of the h ost galaxy (e.g., IMagorrian et al.l II998I : 
IMarconi fc Hunt 2003), indicating the co-evolution be- 
tween SMBHs and galaxies. According to theoretical 
predictions, active galactic nuclei (AGNs) may play a key 
role in regulating both mass accretion and star forma- 
tion by their feedback through powerful outflows and/or 
radiation to the surrounding gas (e.g., IDi Matteo et al.l 
|2005[) . Broadband X-ray observations of AGNs are par- 
ticularly important to study the physics of accretion and 
outflow in the central engine by measuring the spectra 
and time variability of the continuum emission from the 
innermost disk and relativistic jets. They also provide 
unique information on the structure of the tori surround- 
ing the accretion disk, which may be the s ource of mass 
accretion and site of nuclear st arbursts (c.g.. lDavies et all 
120071: llmanishi fc Wadal[200l ). through investigation of 
its reprocessed emission. 

Recent hard X-ray all sky surveys above 10 keV con- 
ducted by Swift/BAT and INTEGRAL are providing us 
with an ideal sample of AGNs in the local universe that 
can be efficiently followed-up with observatories having a 
narrow field-of-view, like Suzaku and XMM-Newton, for 
detailed spectroscopic studies. Suzaku has unique fea- 
tures covering the wide energy band over 0.1-600 keV 
with good energy resolution and large collecting area be- 
low 10 kcV. Making use of these capabilities, we have 
been working on a program of Suzaku follow-up obser- 
vations of Swift / BAT AGNs, mostly focusing on heav- 
ily obscured ones (lEguchi et al.l l200a IWinter et aLll2009l: 



lEguchi et al.l l201lf r A remarkable finding is the dis- 
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covery of "new type" AGNs that exhibit a very small 
amount of scattering gas, which are interpreted to have 
geometrically very thick tori, unlike c lassical type AGNs 
(|Ueda et al.]l2007t lEmichi et al.l l200g). 

Radio galaxies, a class of AGNs exhibiting powerful 
jets, are key objects to understand the origin and effects 
of the AGN feedback to the surroundings. However, the 
fundamental question, how the structure of the nucleus 
is different between AGNs with jets and without jets, 
is not fully resolved yet. This paper is the fourth in a 
series related to Suzaku and Swift/BAT observations of 
AGN, targeting two narrow line radio galaxies (type-2 
radio loud AGNs) 3C 403 and IC 5063, for which no "si- 
multaneous" broad band spectra over the 0.1-200 keV 
band have been reported. All targets studied in our pre- 
vious series of papers are Seyfert galaxies (i.e., AGNs 
without powerful je ts) except for NGC 612 reported in 
lEguchi et al.l (|2011[ ). By adding these new targets, we in- 
crease the sample for comparing the statistical properties 
of radio loud AGNs to radio quiet ones. 

3C 403 {z = 0.059) is a NLRG classified as "X -shapcd" 
or "winged" radio morphology (|Leahv fc Williams T984I) 
, which is explained as the results of propagation 
of je ts in non-axisymmetric atmospheres ([Kraft et al.l 
[2005D The host galaxy Hu bble type is SO or pecu- 
liar. iVasudevan et al.l (|2010( ) obtained the mass of the 
SMBH to be 10®-^^ Mq using the correlation between 
black hole mass and b ulge luminosity, according to the 
procedure described in (jVasudevan fc Fabianll2009l ). The 
bulge luminosity is derived from the total K-band flux 
in the 2MASS point source catalog by subtracting the 
AG N contributi o n est imated from the X-ray luminos- 
ity. iKraft et al.l (120051 ) studied the X-ray properties of 
3C 403 with Chandra, detecting X-ray emission from 
the nucleus, hot interstellar medium (ISM), extended ra- 
dio structures (lobes and wings), and compact regions 
in the jet (knots and hot spots). IC 5063 (z = 0.0110) 
has a highly anisotropic i onizing radiation fie ld with "X" 
or "conical" morphology ([Colina et al.lll99i|) . associated 
with a giant elliptical or SO host galaxy. The black hole 
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TABLE 1 
List of Targets 



Target Name (Swift ID) 



3C 403 (J1952. 4+0237) 



IC 5063 (J2052.0-5704) 



R.A. (J2000)^ 
DecL {J2000)'' 
Redshift^' 


19 52 15.809 


20 52 02.34 


+02 30 24.18 


-57 04 07.6 


0.059 


0.011348 


Suzaku Observation ID 


704011010 


704010010 


Start Time (UT) 


2009-04-08T19:53:00 


2009-04-24T15:04:25 


End Time (UT) 


2009-04-10T02:35:24 


2009-04-25T16:28:19 


Exposure'' (XIS) (ks) 


47.9 


45.2 


Exposure (HXD/PIN) (ks) 


39.3 


57.6 



^ The position and rcdshift for each source are taken from the NASA/IPAC Extragalactic Database. 
Based on the good time interval of XIS-0. 



mass i s estimated to be 10 ''-^^ Mq ijVasudevan et al.l 
l2010f ). iKovama et al.l ()1992D observed this object with 
the Ginga satelhtc and found that the energy spectrum 
in the 3-18 keV band is described by a power-law with 
a photon index of 1.5 and an absorpti on column density 
of 2 X 10^^ cm~2. IVignali et"aLl (|1997[ ) report the results 
of observations with ASCA and ROSAT. From the 0.1- 
10 keV spectrum of IC 5063, they detect two power-law 
components (with photon indices of 1.7±0.2 and 2.2±0.3, 
where the former is absorbed with 2 x 10^'^ cm^^), 

and a narrow iron-Ka line at E'kq! ~ 6.4 keV (with an 
equivalent width of ~200 cV). 

The organization of this paper is as follows. Section 2 
summarizes the observations. Data analysis and results 
are presented in Section 3. We discuss the implica- 
tion of our results in Section 4. The cosmological pa- 
rameters (ifn, ^m, ^ \)= (70 km s'^ Mpc'^ , 0.3, 0.7; 
IKomatsu et al.|[2009() are adopted in calculating the lu- 
minosities. The errors attached to spectral parameters 
correspond to those at 90% confidence limits for a single 
parameter. 

2. OBSERVATIONS 

Suzaku, the fif th Japanese X-ray satellite 
(jMitsuda et al.l I2007D . carries a set of X-ray CCD 
cameras called the X-ray Imaging Spectrometer (XIS) 
and a non-imaging instrument called the Hard X-ray 
Detector (HXD) composed of Si PIN photodiodes 
and Gadolinium Silicon Oxide (GSO) scintillation 
counters. In this paper, we analyze the data of three 
XISs and HXD/PIN, which cover the energy band of 
0.2-10 kcV and 10-60 keV, respectively. The data of 
HXD/GSO covering energies above 50 keV are not 
utilized because the fluxes of our targets are too faint 
to be detected. In the spectral analysis, we also use the 
spectra of Swif t/BAT in the 14-195 keV band averaged 
for 58- months (iBaumgartner et al.ll201ll ). 

The two NLRGs 3C 403 an d IC 5063, cataloged in the 
Swift/B AT 9- months survey (jTueller et al.ll2008D . were 
observed with Suzaku in 2009 April. The basic informa- 
tion of our targets and observation log are summarized 
in Table [TJ Each object was observed at the "HXD nom- 
inal" pointing position for a net exposure of ~50 ks after 
data screening. 

3. ANALYSIS AND RESULTS 

The data reduction was performed according to stan- 
dard procedures from the cleaned event files provided 
by the Suzaku team (processing version 2.3.12.25). Wc 
used FTOOLS (heasoft version 6.8) for extraction of light 



curves and spectra, and XSPEC (version 12. 5. In) for 
spectral fitting. The XIS events of the source were ex- 
tracted from circular regions centered on the source peak 
with a radius of 2' (3C 403) and 2'.5 (IC 5063), in which 
85% and 91% of the total source photons are accumu- 
lated by the X-ray telescopes, respectively. The back- 
ground data were taken from annulus regions centered 
at the averaged optical axis of the X-ray telescopes with 
radii between 2.8'-6.7' for 3C 403 and 1.2'-6.3' for IC 
5063, by excluding regions where the target and other 
bright sources are located. We generated the RMF files 
of the XI S with xisrmfnen, an d the ARF files with xissi- 
marfgen (jlshisaki et al.ll2007D . The "tuned" NXB event 
files provided by the HXD team were utilized to produce 
the background spectra, to which that of the cosmic X- 
ra y background (CXB ) was added based on the formula 
of IGruber et"an JHH). The source flux of HXD/PIN in 
the 15-40 keV band corresponds to ~ 6% (3C 403) and 
~ 24% (IC 5063) of the background (NXB+CXB) rate. 
The signal-to-noise ratio calculated as the source photon 
counts divided by the square root of the total (source 
plus background) counts in this band is 5.3 for 3C 403 
and 26.7 for IC 5063. The statistical error is larger 
than the systematic uncertainty of ~ 0. 34% (Ig) in the 
backg round model for a 40 ks exposure (jFukazawa et al.l 
120091 ). Wc used ac_hxd_pinhxnome5_20080716.rspfor the 
HXD/PIN response. 

3.1. Light Curves 

Figure [T] shows the light curves of 3C 403 and IC 5063 
in the 2-10 keV band combined from XIS-0 and XIS-3 
(upper), those in the 15-40 keV band from HXD/PIN 
(middle), and their hardness ratio between the two 
bands. A bin width of 5760 s, the orbital period of 
Suzaku (including periods of data gap that are excluded 
to calculate the count rate), is chosen to eliminate any 
systcmatics caused by the orbital change of the satellite. 
There is no evident time variability in the light curves 
of 3C 403, while the 2-10 keV fiux of IC 5063 increases 
after ~ 6 x 10^ s from the beginning of the observation. 
The constant flux model is rejected by with reduced 
of 6.34 with 15 degrees of freedom. Due to the limited 
statistics in the PIN data, however, the 15-40 keV light 
curve and hardness ratio of IC 5063 show no significant 
variability based on tests. We examine the spectral 
variability of IC 5063 in Section 13.3.21 

3.2. BAT Spectra 

Before performing spectral fitting to the Suzaku data, 
wc analyze the Swift/BAT spectra in the 14-195 keV 
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Fig. 1. — (top) Light curves of 3C 403 obtained with the XIS in 
the 2-10 keV band (upper) and with the HXD/PIN in the 15-40 
keV band (middle), and the hardness ratio between the two bands 
(lower). The bin size is 5760 s. The reduced with the degrees 
of freedom for the constant flux hypothesis is shown at the upper 
right corner in each panel, {bottom) The same as IC 5063. 



band to constrain the high energy cutoff. The de- 
tection significance of 3C 403 and IC 5063 in the 58- 
months data is 9.6 and 26.8, respectively (Baumgart- 
ner et al. 2011). A simple power law fit with oc 
is statistically acceptable, yielding a photon index of 
r 1.7 ± 0.3 with xVdof = 8.5/6 for 3C 403, and 
r = 1.9 ± 0.1 with xVdof = 3.7/6 for IC 5063. It 
is known, however, that the hard X-ray continuum of 
AGNs is generally described by a power law with an 
exponential cutoff, cx i?~'"cxp(— _B/_Bcut), where i?cut 
is a cutoff energy. S ince the cutoff energy is typi- 
cally > 100 kcV (e.g.. iDeluit fc Courvoisieil 120031 ). the 
Suzaku data alone (XIS and HXD/PIN) cannot deter- 
mine it. To fit the BAT spectr a, we adopt pexrav model 
(|Magdziarz fc Zdziarskil995() in XSPEC, which contains 
a Compton reflection component from optically thick 
cold matter. The strength of reflection is represented 
with R = Jl/27r, where i7 is the solid angle of the reflec- 
tor viewed from the irradiating source. In this analysis, 
we assume two extreme cases for the reflection strength, 
R = and 2. The inclination angle in the pexrav model 
is fixed at 60° as a representative value throughout our 



TABLE 2 

Cutoff Energies (Bcut) Estimated by the BAT spectra 





Target Name 


3C 403 


IC 5063 


R = 


Scut 

xVdof 


540(> 120) 
8.9/6 


220lgo° 
3.3/6 


R = 2 


Ecut 

xVdof 


570(> 140) 
8.6/6 


12.9/6 



Note. — The photon index is fixed at 1.7 for both targets. 



paper except in Section 13.41 To avoid strong coupling 
between i^cut and F, we adopt F = 1.7 for both targets, 
which is an averaged value of 13 near by radio galaxies ob - 
taincd from previous X-ray studies (jCrandi et al.|[2006l) . 
As summarized in Table [21 we find that Ecut > 120 keV 
for 3C 403 and Scut = 130-470 keV for IC 5063. In the 
following spectral analysis, we fix the cutoff energy at 200 
keV for both 3C 403 and IC 5063, since changing these 
values within the above uncertainties does not affect our 
conclusions. 

3.3. Spectral Analysis with Analytic Models 

For spectral analysis, we use the Suzaku spectra of the 
2 XIS-FIs in the 1-11 keV band, and the XIS-BI in the 
0.5-8 keV band, and HXD/PIN in the 16-40 keV band, 
where the highest signal-to-noise ratios are achieved. The 
XIS data in the 1.7-1.9 keV band are discarded because 
of calibration uncertainties associated with the instru- 
mental Si-K edge. The relative flux normalization of 
the XIS-BI to the XIS-FIs is set free, while that of the 
HXD/PIN to the XIS-FIs is fixed at 1.18 based on th e 
calibration using the Crab Nebula (IMaeda et al.ll2008D . 
The Galactic absorption, N^'^^, is included in the spec- 
tral model, which is fixed a t the v alue derived from the 
H I map of iKalberla eFall (|2005l) {N§''^= 1.22 x lO^i 
cm-2for 3C 403 and 6.10 x 10^° cm-^for IC 5063). Solar 
abundances by lAnders fc Grevess^ ()1989[ ) are assumed 
throughout our analysis. The upper panels of Figure [5] 
plot the Suzaku spectra of the two sources folded with 
the detector responses, together with the BAT spectra 
in units of photons cm~^ s~^. 

Wc apply the three analytic models defined in 
lEguchi et all ()2009[ ). Models A, B, and C (see below), 
to which we add a spectral component from an optically 
thin plasma with Solar abundances (apec) because it 
improves the fit significantly in both targets. In the case 
of 3C 403, two power law components from the radio 
lobes/ wings and from the compact structures in the jets 
are also included. These ~10" scale structures are spa - 
tially resolved in the Chandra data (jKraft et al.l 12005} ); 
their fluxes in the 0.25-10 keV band are estimated to be 
(4.4 ±0.3) X 10-14 erg cm'^ s'^ and (3.5 ±0.3) x 10^" 
erg cm^^ s^^, corresponding to 0.31% and 0.27% of the 
total flux from the unresolved nucleus corrected for ab- 
sorption, respectively. 

Model A consists of a direct component from the nu- 
cleus (a cutoff power law absorbed by cold matter), a 
scattered component (a cutoff power law without absorp- 
tion), and an iron-K emission line (a Gaussian). The 
high energy cutoff is fixed at 200 keV based on the BAT 
only results (see above). The scattered component is 
assumed to have the same slope as the direct one with 
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a fraction of /scat- In Model B, a Conipton reflection 
component of the direct continuum from optically thick 
matter (pexrav model), absorbed with a column den- 
sity of A^f», is further added to Model A. This emission 
is expected from the inner wall of the torus and/or the 
accretion disk irradiated by the central source. The ab- 
sorption for the reflection component is set to be indepen- 
dent of that of the transmitted one, because the emission 
region is different and the line-of-sight column density 
may not always be the same, depending on the viewing 
geometry. Model C has the same emission components 
as Model B including an absorbed reflection component, 
but we consider two layers of absorber with different 
hydrogen column densities for the transmitted compo- 
nent, corresponding to the "partial covering" case where 
the absorbers in the line of sight consist of gas blobs 
smaller than the size of t he X-ray emitting region (e.g., 
lElitzur fc Shlosmanl[200fih . 

We flrst perform simultaneous fit only to the Suzaku 
XIS and HXD/PIN spectra to select the most appropri- 
ate model for each object among the three. To avoid 
complex effects of time variability, the BAT spectra av- 
eraged for 58 months are not utilized in this stage. We 
start from the simplest model (Model A), and adopt a 
more complex model (Models B or C) only if we find a 
statistically significant improvement of the fit or a phys- 
ically more reasonable solution. After selecting the best 
model describing the data in this way, we finally include 
the BAT spectra in the joint fit to better constrain the 
continuum up to 200 keV. We allow only the normaliza- 
tion of the direct component to vary while keeping the 
continuum shape (i.e., spectral slope and cutoff energy) 
the same between the Suzaku and Swift/B AT epoch. The 
normalizations of the refiection and scattered component 
are also linked between Suzaku and Swift/HAT, consid- 
ering that the time scale of their variability should be 
larger than >years if the emission region size is typical 
(>1 pc) of the scale size of tori. Unlike our prev ious 
papers (|Ueda et all [20071: lEguchi et al1[200l [20Tll ). we 
define the refiection strength and scattered fraction rel- 
ative to the BAT flux, not to the Suzaku flux, since the 
BAT spectra integrated over 58 months should be a good 
indicator of the averaged flux level responsible for these 
reprocessed emission. 

3.3.1. 3C 403 

From the results by iKraft et al.l (|2005| ) obtained with 
Chandra, we are able to estimate the X-ray fluxes from 
the radio lobes/wings and compact features in the jets 
(knots and hot spots) that should contribute to the 
Suzaku spectra, where these jet-related structures and 
the nucleus are not spatially resolved. Since neither of 
these components are expected to be variable on time 
scale of years, we fix their fluxes at the best-fit Chan- 
dra values; we confirm that the uncertainties do not af- 
fect our results as their contributions in the total flux 
are minor. The component from the radio lobes/ wings 
and jets are represented as a power law with a photon 
index of 1.8 and 2.0 with a normalization at 1 keV of 
9.8 X 10~^^ and 9.6 x 10~^^erg cm~^s~^keV~^, respec- 
tively, both are always added to the above three mod- 
els representing the nucleus emission. Since photon in- 
dex of the power law in the lobes/ wings component is 
poorly constrained from the Chandra data, we fix its 
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Fig. 2. — (top) The folded spectra of 3C 403 obtained with the 
XIS-FIs in the 1-11 kcV band (black), the XIS-BI in the 0.5-8 keV 
band (red, open circles), the HXD/pIN in the 16-40 keV band 
(blue, open squares), and the Swift/HAT in the 14-195 kcV band 
(magenta, open triangles). The best-fit model (B in Table [Sj is 
plotted by the solid lines, and the residuals in units of x a-re shown 
in the lower panels, (bottom) The same as IC 5063. 



slope to that of the radio emission iDennett-Thorpe et al.l 
l|2002f ). assuming that the radio and X-ray emission cor- 
respond to the synchrotron and inverse Comptonization 
from the same electrons, respectively. To model optically 
thin thermal emissi on from the hot in terstellar medium 
(ISM) detected bv IKraft etall ()2005f ). we also add an 
apec model whose temperature and emission measure 
are free parameters. The elemental abundances are fixed 
at the Solar values. 

We adopt Model B as the most appropriate model 
for the nucleus emission of 3C 403. Model A yields 
a very small photon index of the transmitted compo- 
nent, '--jI.O, which seems rather unphysical, with x^/dof 
= 82/73. Reasonable results arc obtained with Model B, 
which gives a photon index of ~1.5 with x^/dof = 74/71. 
Model C does not improve the fit significantly in terms 
of the value (x^/doi = 74/69). The best fit model 
(Model B 4- two power laws -I- apec) to the joint fit in- 
cluding the Swift/'BAT spectrum is overplotted in the up- 
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Fig. 3. — (top) Best-fit analytic model for 3C 403 in units of 
EF(E) {F{E) is tlie energy flux at the energy E). The solid 
(black), dashed (red), dotted (blue), and dot-dashed (magenta) 
curves correspond to the total, direct component, reflection com- 
ponent, and iron-K emission line. The dashed (purple) curve rep- 
resents the emission from an optically thin thermal plasma, and 
dot-dot-dot-dashed (orange) curves show the emission from the 
jets and radio lobes/wings (see the text). The dot-dot-dot-dashed 
(cyan) curve represents the scattered component, (bottom) The 
same as IC 5063. 



the black hole, fully consistent with a torus origin. The 
equivalent width of the iron-K line with respect to the 
reflection continuum is 1.13±0.24 keV, which is physi- 
cally self-consistent acc ording to theoretical calculations 
fe.g.. llkeda et al.|[2009l ). The smaller column density for 
the reflected component {N^^ 4 x 10^^ cm~^) than 
that for the transmitted one (iVn^ 6 x 10^"^ cm"^) sug- 
gests that we are looking at the source with a somewhat 
face-on geometry so that a significant fraction of the re- 
flected light from the inner wall of the torus can reach us 
without strong obscuration. 

While we detect a very strong reflection within the 
Suzaku spectra corresponding to i? « 2, which dominates 
the HXD/PIN flux, we find a lower value of i? « 0.6 when 
normalized to the flux of the transmitted component de- 
termined by Swift/BAT. This can be explained by time 
variability since the flux during the Suzaku observation 
was about 3 times smaller than the averaged flux over the 
58 months of the BAT observations. In fact, the differ- 
ence between the unabsorbcd nuclear flux in the 0.25-10 
keV band was 1.3 x 10"^^ erg cm~^ s~^ at the Chandra 
observation in 2002, while that measured with Suzaku is 
4.0 X 10~^^ erg cm~^ s~^. 

The best fit temperature of the apec component is 
kT = 0.29 ±0.03 k eV, consistent with the Chandra result 
(|Kraft et al.ll2(105h . kT = 0.3 ± 0.04 keV. Wc note, how- 
ever, that its flux as determined with Suzaku, 2.8 x 10^^'^ 
erg cm~^ s~^ in the 0.25-10 keV band, is about 5 times 
higher than that estimated from the Chandra data, even 
taking into account the contribution from the hot ISM in 
an outer region at radii larger tha n 7". 5 up to 60" from 
the nucleus (see lKraft et al.l ()2OO50 Section 3.1), which is 
contained in the Suzaku extraction region. We find that 
this discrepancy is mainly produced due to the spectral 
coupling with the unabsorbed soft component; its slope 
is assumed to be the same as for the direct component in 
our spectra l model, F = 1.5, w hich is smaller than that 
adopted bv iKraft et all ()2005D . F = 2.0. Thus, the dif- 
ference may be regarded as a systematic uncertainty in 
the absolute flux from the hot ISM, whose core emission 
is difficult to spatially resolve even with Chandra from 
the AGN component. We confirm that our result on the 
amount of the scattering gas presented in Section 13.3.11 
is not significantly affected (within a factor of 1.5) even 
when we fix the flux of the apec component at the Chan- 
dra value in the spectral analysis. 



per panel of Figure [2] with residuals in the lower panel. 
Figure [3] (a) plots its EF{E) form, where F{E) is the 
energy flux at the energy E. 

The best-fit parameters are summarized in Table [3l 
The moderately flat photon index (F 1.53j;o o2): 
though within, the typical range for radio galaxies 
(IGrandi et al.l 120061) . may be affected by simplification 
of the analytic model where only a single absorber is 
assumed for the reflection component. The scattered 
fraction /scat ~ 0.4% is rather small, corre sponding to 
that of a "new type" AGN in lUeda et al.l (|2007) if it 
corresponds to the opening angle of the torus. This is- 
sue will be discussed in Section 4. The iron-K emission 
line is detected at « 6.4 keV and is narrow (with a la 
width of <98 eV at 90% confidence level). This indicates 
that it mainly comes from cold and distant matter from 



3.3.2. IC 5063 

Model B with thin thermal emission (apec) is adopted 
as the best model for IC 5063 from the time averaged 
Suzaku spectra. The thermal emission is required be- 
cause it significantly reduces in all 3 models (A, B, 
and C) and its addition produces acceptable fits. Mod- 
els A, B, and C give xVdof = 185.8/177, 175/175 and 
175/173, respectively. Thus, the improvement of the fit 
by adding a reflection component (Model B) is found to 
be significant with an F-test probability of 0.007, whereas 
the introduction of an additional absorber for the trans- 
mitted component (Model C) is not significant in intro- 
ducing . The results of simultaneous analysis of the 
Suzaku and BAT spectra with Model B + apec are sum- 
marized in Table [3l and the best-fit model is plotted in 
Figures H and H 
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TABLE 3 

Best-Fit Parameters of Suzaku and Swift/BAT spectra with Analytic Model 





Target 


3C 403 


IC 5063 




Best-fit model 


B -|- apec -|- two power laws 


B + apec 


(1) 


^Gal(iQ22 


0.122 


0.061 


(2) 


A^hCIO^^ cm 


-r.-! -1-6 

6II5 


25.0 ± 0.7 


(3) 


NormBAT 


3.2 ± 0.7 


1.06 ± 0.08 


(4) 


r 


i_n 

J^-a'J_0.02 


1.72 ± 0.02 


(5) 


/scat (%) 


0.4 ±0.3 


0.9 ±0.1 


(6) 


Econ (keV) 




ft QQ+0.02 

"•■^^-0.08 


(7) 


£^wid (eV) 


< 98 


< 78 


(8) 


EW (eV) 


455 ± 98 


148 ± 20 


(9) 


EW^*' (keV) 


1.13 ±0.24 


1.52 ± 0.21 


(10) 


Ag=fl(1022 cm-2) 




3 1+3.4 


(11) 


R 


0.65 ±6.17 


0.65 ± 0.13 


(12) 


kT (keV) 


0.29 ± 0.03 




(13) 


n'^V (cm"'') 


(2.7 ±0.7) X 10^5 


(1.0 ±0.2) X 10*53 


(14) 


F2-10 (erg cm~2 8"^) 


7.14 X 10-13 


8.03 X 10-12 


(15) 


(erg cm~2 s^^) 


7.96 X 10-12 


3.76 X 10-11 


(16) 


L2-10 (erg s^i) 


1.76 X 10"^ 


5.61 X 10*2 






81.7/78 


182.6/182 



Note. — Errors correspond to 90% confidence level for a single parameter. Both /scat and -R are normalized to the direct component 
flux as measured with BAT. 

(1) The hydrogen column density of Galactic absorption bv lKalberla et aP l|2005l 'l. 

(2) The linc-of-sight hydrogen column density for the direct component. 

(3) Normalization ratio of the direct component between the BAT and Suzaku spectra. 

(4) The power-law photon index of the direct component. 

(5) The fraction of the scattered component relative to the intrinsic power law. 

(6) The center energy of the iron-K« emission line at the rest frame. 

(7) The observed width of iron-Ka line. This is fixed at 1 eV when we determine the continuum parameters. 

(8) The observed equivalent width of the iron-Ka line. 

(9) The observed equivalent width of the iron-Ka line with respect to the reflection component. 

(10) The line-of-sight hydrogen column density for the reflection component. 

(11) The relative strength of the rcficction component to the direct one, defined as /? = f2/27r, where Q is the solid angle of the rcfiector. 

(12) The temperature of the apec component. 

(13) The emission measure of the apec component. 

(14) The observed Suzaku fiux in the 2-10 keV band. 

(15) The observed Suzaku fiux in the 10-50 keV band. 

(16) The 2-10 kcV intrinsic luminosity obtained with Suzaku, corrected for the absorption. 



The best fit photon index of F w 1.7 from IC 5063 is 
typical of radio galaxies. The scattered fraction /scat ~ 
0.9% is larger than that o f 3C 403. It is cons istent with 
the ASCA result of ~1% (jVignah et al.lflQQTl) . The nar- 
row iron-K line detected at 6.4 keV is consistent with 
being emitted from the torus. Its equivalent width with 
respect to the reflected continuum (1.5 ± 0.2 kcV) is 
within a range of theoretical expectation, thus support- 
ing the validity of our spectral model. Similar to the 
case of 3C 403, the reflection component is less absorbed 
{N^^ « 3 X 10^^ cm-^) compared with the transmit- 
ted component (iVH~ 3 x 10^'^ cm-^), and disfavoring 
the likely edge-on geometry. The reflection strength of 
R « 0.6 relative to the averaged BAT flux is close to the 
value obtained from 3C 403. These results suggest that 
the torus geometry of the two AGNs could be similar to 
each other. The obtained temperature and luminosity of 
the apec component are in the t ypical range of hot ISM 
detected from elliptical galaxies (jlVIatsushita 
and hence supports this interpretation of its origin, al- 
though we cannot rule out the possibility that a part of 
this emission may come from photoioinized plasma irra- 
diated by the AGN, which should be regarded as an ac- 
companying component of the scattered emission. Nev- 
ertheless, since the amplitude of the scattered fraction is 
mainly determined from the spectrum in 1-2 keV band, 



where no strong emission lines are present, its intensity 
is not strongly affected by the uncertainty in the nature 
of the line emitting component. 

Since IC 5063 exhibits significant time variability dur- 
ing the Suzaku observation, we check time variability of 
the spectra by dividing the time region into two states 
with different flux levels, epoch I (for 57.6 ks from the be- 
ginning) and epoch II (after that). We adopt the Model 
B + apec model for each spectrum. Assuming that the 
reflection component and thin thermal emission were not 
variable on this short time scale, since they most proba- 
bly originate from the torus and host galaxy, respectively, 
we fix the spectral parameters of these components (and 
its absorption) at the best-fit values determined by the 
time-averaged Suzaku (and Swift/BAT) spectra. We find 
that the shape of the transmitted component did not 
show variability within statistical errors, T = 1.7 ± 0.1 
andA^H = (26±2)xl022 cm-^in epoch I and T = 1.8±0.1 
and A^H = (25 ± 2) x 10^^ cm-^in epoch II, and only the 
normalization changed between the two epochs. Thus, 
averaging the Suzaku spectra over the whole epoch is 
justified. 

3.4. Application of Torus Model 

We also fit the spec t ra gen erated by the Monte-Carlo 
model of llkeda et all (|2009D . where reflection compo- 
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TABLE 4 

Best-Fit Parameters of Suzaku and Swift/BAT spectra with Torus Model 







3C 403 


IC 5063 


(1) 


Model 


Solar + apec" + two power laws 


Solar + apec'' 


(2) 


AfGal(io22 cm-2) 


0.122*' 


0.061'' 


(3) 


N^^ (1022 cm-2) 


67 ±25 


25if 


(4) 


^oa (degrees) 


50'' 


50'' 


(5) 


fine (degrees) 






(6) 


NorniBAT 


3.2 ± 1.2 




(7) 


r 


, r-c+0.32 




(8) 


/scat,o(%) 


0.3(< 0.7) 


0.5 ±0.1 


(9) 


Eccn (keV) 


6.4*' 


6.4'' 


(10) 




1 2+°-'' 


-1 n + 0.2 
^■^-0.3 






88.4/80 


182.0/187 



Note. — Errors are 90% confidence level for a single parameter, /scat is normalized to the direct component flux as measured with 
BAT. 

(1) The model used in the fit. "Solar" means the torus mo del with Solar abundan ces. 

(2) The hydrogen column density of Galactic absorption bv lKalberla et aTl II2005I ). 

(3) The hydrogen column density of the torus along the equatorial direction. 

(4) The half opening angle of the torus. 

(5) The inclination angle of the torus. The lower limit is restricted to be 51° by the condition 8^^^ > 0oa + 1. 

(6) Normalization ratio of the direct component between the BAT and Suzaku spectra. 

(7) The power-law photon index. It is restricted in the range of 1.5-2.5 in the table model. 

(8) The fraction of the scattered component relative to the intrinsic power law when the half opening angle of the torus is 45° . 

(9) The center energy of the iron-K emission line at the rest frame of the source redshift. 

(10) The relative strength of the iron-K emission line to that predicted by the torus model. 
The parameters of the apec model arc fixed at those in Table 3. 

'■ The parameters are fixed at this value. 



ncnts from a three-dimensional, axially symmetric torus 
around the nucleus emitting a cutoff power law contin- 
uum is computed. The torus has four parameters, the 
half opening angle (^oa): inclination {O-mc), total column 
density in the equatorial plane (N^), and the ratio be- 
tween the outer and inner rad ii, which is f i xed to be 100 
(for definition, see Figure 2 of llkeda et all (|2009[ )). 

We analyze the data in the same manner as described 
in lEguchi et al.l (j2011l ) , to which we refer the reader for 
details of the spectral formula in the XSPEC terminol- 
ogy. The unabsorbed and absorbed reflection continua 
and fluorescence iron-K line from the torus are included 
in the model as separate components. The line-of-sight 
column density of the direct component is automati- 
cally determined from the four geometrical parameters 
of the torus; for absorbed AGNs 9oa < dine- The ob- 
served scattered component is represented as /scat = 
/scat,o(l -cos6'oa)/(l -cos45°), where /scat.o reflects the 
averaged column density of the scattering gas. The re- 
flection component from the accretion disk is also consid- 
ered for the direct component by assuming a solid angle 
of f2 = 27r with the same inclination as for the torus. 
Solar abundances are assumed. We fit the Suzaku (XIS 
and HXD/PIN) and Swift/BAT spectra simultaneously 
by limiting the energy band to 0.5-100 keV to which 
the Ikcda model is applicable. Although the model uses 
a cutoff energy of 360 keV instead of 200 keV, the ef- 
fect is negligible as we only use energies below 100 keV. 
Similar to the analysis in Section [3.31 time variability of 
the flux of the direct component between the Suzaku and 
Swift/BAT epochs is taken into account with an assump- 
tion that the power law slope is constant. The apec com- 
ponent and the two power law components originating 
from the radio lobes/wings and resolved jet components 
(for 3C 403) are included in the spectral model, whose 
parameters are fixed at the values found in Section [3. 3. II 



We find that this "torus model" also provides a good 
fit to the observed spectra of 3C 403 and IC 5063 with 
;^2/dof=88/80 and 182/187, respectively. The torus 
opening angle is not well determined in the model due to 
the limited photon statistics and we fixed it at 6'oa = 50° 
in our analysis consistent with the optical observations of 
the ionized cone in IC 5063 (- 50°; [Colina ct al. 1991). 
The free parameters of the spectral model are 6'i„c, Af^, 
/scat, J the normalization ratio between the iron-K emis- 
sion line and reflected continuum (epe), the photon index 
(F), and two normalizations of the cutoff power law for 
the Suzaku and Swift/BAT data. The best-fit parameters 
are summarized in TablcUl Wc find /scat,o = 0.3(< 0.7)% 
(1.0l^J%) for 3C 403 and /^cat.o = 0.5±0.1% (0.6l°;^%) 
for IC 5063 when normalized to the BAT (or Suzaku) 
fiux. These results on the scattering gas will be discussed 
in Section 4. In both targets, the best-fit photon index 
becomes slightly larger than in the analytical model, al- 
though the difference is within the statistical errors. The 
normalization ratio of the direct component between the 
BAT and Suzaku spectra is consistent with that obtained 
from the analytic model. As discussed in the previous 
subsections, the data indeed favor an inclination angle 
that is only slightly larger than the torus opening angle, 
particularly for 3C 403. 

A si milar numerical torus model bv lMurphv &: YaqoobI 
l|200g ) (MYTORUS in XSPEC) also gives a good fit for 
both targets, where a slightly different torus geometry is 
assumed with a fixed opening angle of 60° . By adopting 
it instead of the Ikeda model in the same way, we obtain 
F = 1.7l°;i, 0i„c = 75t}° deg, and N^^ = O.StH x 10^^ 
cm-2for 3C 403 with xVdof=85.5/79, and F = 1.9 ± 
0.1, dine = 63l^ deg, and N^^ = 0.6t°| x 10^^ cm'^for 
IC 5063 with xVdof=182.3/183. These parameters are 
consistent with those obtained by the Ikeda model within 
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the statistical errors, except for of IC 5063. 

4. SUMMARY AND DISCUSSION 

With Suzaku, we have for the first time obtained si- 
muhaneous broad band X-ray spectra from two NLRGs, 
3C 403 and IC 5063, with CCD energy resolution below 
10 keV. Combined with the Swift/BAT spectra averaged 
for 58 months, we are able to best constrain the spec- 
tral properties over the 0.5-200 keV band. The bolomet- 
ric luminosities of 3C 403 and IC 5063 are estimated as 
10^^'^ erg s~^ and lO^*^'^ erg s~^ from the infrared and 
X-ray luminosities (^V asudevan et al.1l2010D . correspond- 
ing to the Eddington fractions of appro ximately 0.06, us- 
ing th e black hole masses estimated bv lVasudevan et"an 
(j2010t ). These Eddington fra ctions are in the range ob- 
served for radio galaxies (e.g.. lEvans et al.ll2006[) and for 
Seyfert galaxies ([Winter et al.l2009() in the local universe. 
Our targets thus significantly increase the NLRG sample 
with similar Eddington fractions studied with Suzaku and 
Swift/BAT in addition to NGC 6 12, which has Lboi/-^Ed d 
= 0.02 and iboi = lO'^"''^ erg a'^ (jVasudevan et al.|[20l^ . 

While Suzaku provide us with unique spectroscopic 
data, the integrated spectra from radio galaxies in- 
evitably contain the contribution from various compo- 
nents besides the nucleus, such as the hot ISM, the radio 
lobes/wings, and the kpc-scale jets, due to the limited 
spatial resolution. The thermal emission is important 
only in the softest energy band in the Suzaku data. In 
the case of 3C 403, we are able to separate those from the 
radio lobes/wings and k pc-jets, utilizing the Chandra re- 
sults (jKraft et al.ll2005[) . We find that the relative fluxes 
of these components are 1.0% and 0.7% of the unab- 
sorbed nucleus emission at 5 keV, respectively, and thus 
make only minor contribution in the hard X-ray spectra. 
Although we cannot do the same analysis for IC 5063 and 
NGC 612 due to the lack of Chandra data, we estimate 
that the inverse Compton component originating from 
the radio lobes is 0.1% and 1% of the nucleus emission 
at 5 keV, respec tively, u sing the observed radio spectra 
from the lobes ([Jones fc McAdam 1992) and a typical 
X-ray (inverse Compton) to radio (synchrotron) flux ra- 
tio of the lobe emission ([TsobQi2002il . One should note, 
however, that emission from the "nucleus" may also con- 
tain components fro m the pc-scale jets , which are seen 
in th e VLBI images ([Tarchi et al.ll2007l: lOosterloo et al.l 
20001) but cannot be r esolved with Chandra. In fact, 
Evans et al.) ([2OO60 and iHardcastle et al.l (I2009D suggest 
that the unabsorbed soft X-ray components of radio 
galaxies is mainly attributable to the unresolved jets, 
based on the correlation between the soft X-ray and radio 
core luminosities at 5 GHz. In this context, the "scat- 
tered fraction" we discuss in our paper is only an upper 
limit. 

For 3C 403, we can compare our results of the 
scattered fraction and absorption column density of 
the transmitted component obtaine d through X-ray 
spectr oscopy with those estimated by iMarchesini et al.l 
(l2005l) fro m the near-infrare d (NIR ) and optical imag- 
ing data. IMarchesini et al.) ([2005D calculate the scat- 
tered fraction and extinction from the "observed" ra- 
tios of the optical (continuum or line) luminosity to 
the NIR one, assuming the typical "intrinsic" luminos- 
ity ratios obtained for radio loud AGNs. Our result, 
/scat = (0.4 ± 0.3)%, is smaller than their estimate. 



0.88l[j g%, though within the errors. This possible dis- 
crepancy may be explained if there is another unabsorbed 
component in the NIR/optical lights in addition to the 
scattered one, such as thermal emission from hot dust in 
the torus, and/or if there was time variability in the di- 
rect NIR/optical fluxes, which were measured from short 
term observations. The column density derived from 
the X-ray data is '^lO times larger than that converted 
from the estimated extinction. Ay ~ 21.9 mag, assum- 
ing the Ga lactic standard relation Ay = 5 x 10~^^iVH 
([Bohlin et a l. 1978). This suggests much smaller dust- 
to-gas ratio in the torus th an Galactic gas, as r eported 
for many other AGNs (e.g.. iMaiolino et al.lf200l . 

Applying analytic spectral models that are commonly 
used to represent the X-ray nuclear spectra of AGNs (i.e., 
a cutoff power law and reflection component), we con- 
strain the reflection strength in terms of the solid angle 
of the reflector to be i? = fi/27r « 0.6 for both 3C 403 and 
IC 5063. A similar value (i? « 0.5 when normalized to 
the BA T flux) is obtained for NGC 612 bv lEguchi eta[\ 
([20 lit) , while a weaker reflection (i? < 0.4 at 90% confl- 
dence limit) is found from another NLRG, 3C 33. These 
results suggest that a major fraction of NLRGs show only 
a moderately strong reflection compo nent in their X -ray 
spectra, like many Seyfert 2 galaxies (lDadinall2007D . im- 
plying that their tori have a sufficiently large solid angle 
to produce the reprocessed emission. Thus, there is no 
evidence that the torus structure is systematically differ- 
ent between radio galaxies and Seyfert galaxies, although 
studies of a larger NLRG sample are required to inves- 
tigate this issue. The reason for the weak reflection in 
3C 33 is unclear but may be connected to its high X-ray 
luminosity (L ^ 4 x lO'^^ erg s~^ in the 2-100 keV band) 
compared with the other NLRGs. 

The application of the numerical torus model where 
a realistic geometry is assumed gives useful insights 
into the torus str ucture, as demonstrated in our pre- 
vious studies (e.g llkeda et al.l 120091: lAwaki et all 120091: 
lEguchi et al1l201ll) . For our targets, we cannot constrain 
the torus opening angle due to the limited photon statis- 
tics; in 3C 403 since values of 6'oa = 30°, 50°, or 70° yields 
a very similar value, and thus we cannot determine if 
the torus of 3C 403 is close to that of " new type" AGNs or 
to "classical type" ([Ueda et al.ll2007t lEguchi et ani2009[) 
from our data alone. As for IC 5063 we have a reasonable 
constraint on the torus geometry from the measurement 
of the ionization cone, 6'oa = 50° . Together with the fact 
that NGC 612 has 6'oa > 58° ([Eguchi et al.ll20Tl . we 
do not find any evidence for the presence of "new type" 
AGNs in NLRGs studied with Suzaku so far, although 
wc cannot rule out its possibility for 3C 403. We need 
a large sample with good quality broad band spectra to 
verify whether or not the distribution of population is bi- 
mo dal with "new type" and " classical type" as proposed 
by lEguchi et al.l ([2009() , both for radio quiet and loud 
AGNs. 

A notable feature seen in 3C 403, IC 5063, and 
NGC 612, is the small amount of scattering gas, parame- 
terized as /scat, 1 reflecting the averaged column density 
of the scattering gas. We find /scat.o < 0.7% in these 
radio galaxies, which is systematically smaller than that 
obtained for Seyfert galaxies fitted with the same model, 
> 2.0% for NGC 3081 (lEguchi et al.ll20H and 1.3t^ j% 
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for SWIFT J0255. 2-0011 (|Eguchil [20lTI ). in which the 
scattered component is normahzed to the flux as mea- 
sured with BAT. As mentioned above, one must take 
this value only as an upper limit considering the possi- 
ble contribution from the unresolved jets to the soft X- 
ray component. The /scat,o value depends on the torus 
opening angle assumed, but the above conclusion is un- 
changed unless Oos. ^ 30° for 3C 403. Our result suggests 
that there is intrinsically little "normal" gas around the 
nucleus in radio galaxies. This result may be explained 
by the jet activity which can expel such gas. 

We note that the obtained reflection strengths and 
scattered fractions in 3C 403 and IC 5063 arc normalized 
to the hard X-ray flux over 58 months obtained through 
the Swift/BAT monitoring, which should give a good 
measure of the "averaged" activity of each source. In 
fact, without the Swift/BAT data, we obtain an appar- 
ently very strong reflection strength R ^ 2 from 3C 403. 
We demonstrate that such an unphysically large R value 
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